arXiv: 1505.080lOvl [math.CO] 29 May 2015 


Interlacing Families IV: Bipartite Ramanujan Graphs of All Sizes * 


Adam W. Marcus Daniel A. Spielman Nikhil Srivastava 
Crisply Yale University UC Berkeley 

Yale University 


June 1, 2015 


Abstract 

We prove that there exist bipartite Ramanujan graphs of every degree and every number 
of vertices. The proof is based on analyzing the expected characteristic polynomial of a union 
of random perfect matchings, and involves three ingredients: (1) a formula for the expected 
characteristic polynomial of the sum of a regular graph with a random permutation of another 
regular graph, (2) a proof that this expected polynomial is real rooted and that the family of 
polynomials considered in this sum is an interlacing family, and (3) strong bounds on the roots 
of the expected characteristic polynomial of a union of random perfect matchings, established 
using the framework of finite free convolutions introduced recently in IMSS15al . 
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1 Introduction 


Ramanujan graphs are undirected regular graphs whose nontrivial adjacency matrix eigenvalues 
are asymptotically as small as possible; in other words, they are the optimal spectral expander 
graphs. In this paper, we prove the existence of bipartite Ramanujan graphs of every degree and 
every size. We do this by showing that a random m—regular bipartite graph, obtained as a union 
of m random perfect matchings across a bipartition of an even number of vertices, is Ramanujan 
with nonzero probability. Specifically, we prove that the expected characteristic polynomial of such 
a random graph has roots concentrated in the appropriate range, and use the method of interlacing 
families introduced in [MSSI5b] to deduce that there must be an actual graph whose eigenvalues 
are no worse than the roots of this polynomial. Infinite families of bipartite Ramanujan graphs 
were shown in that paper to exist for every degree m > 3 , but it was not known whether they 
exist for every number of vertices. 

The main conceptual and technical contributions of this work and the companion paper [MSS 15a] 
are the following. First, we identify a new class of real-rooted expected characteristic polynomials 
related to random graphs, and develop new tools for establishing their interlacing properties and 
analyzing the locations of their roots. These methods are different from those used to study the 
mixed characteristic polynomials of [MSS 15c] . and the bounds we obtain are strictly stronger than 
those produced by the original “barrier method” argument introduced in BSS121 (which is off 
by a factor of two in this setting). Notably, the expected characteristic polynomials we consider 
are computable in polynomial time, unlike most other known expected characteristic polynomials. 
Second, in contrast to previous work, we derive the Ramanujan bound from completely generic 
considerations involving random orthogonal matrices, in particular making no use of results from 
algebraic graph theory or number theory. 


1.1 Summary of Results 

Recall that the adjacency matrix A of an m—regular graph on d vertice^j] has largest eigenvalue 
Ai(A) = m, and smallest eigenvalue \d{A) = —m when the graph is bipartite. Following Friedman 
[FriOBj . we will refer to these as the trivial eigenvalues of A, and we will call a graph Ramanujan if all 
of its non-trivial eigenvalues have absolute value at most 2 \Jm — 1. Such graphs are asymptotically 
best possible in the sense that a theorem of Alon and Boppana [Nil91] tells us that for every e > 0, 
every infinite sequence of m —regular graphs must contain a graph with a non-trivial eigenvalue of 
absolute value at least 2 \Jm — 1 — e. 

Our main theorem is that a union of m random perfect matchings across a bipartition of 2 d 
vertices is Ramanujan with nonzero probability. 

Theorem 1.1. Let Pi,... , P m be independent uniformly random dxd permutation matrices, m > 3. 
Then, with nonzero probability the nontrivial eigenvalues of 


-4 = £ 



Pi 

0 


are all less than 2 y/m — 1 in absolute value. 

1 In order to be consistent with our companion paper IMSS15aj . we will, unconventionally, use m to denote the 
degree of a graph and d to denote its number of vertices. 
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We also prove the following non-bipartite version of this theorem, regarding a union of m random 
perfect matchings on d vertices (not bipartite), with d even. 

Theorem 1.2. Let P\,.... P rn be independent uniformly random dx d permutation matrices, d 
even, m > 3. Let M be the adjacency matrix of any fixed perfect matching on d vertices. Then 
with nonzero probability: 



Since we only prove nonzero bounds on the probabilities, the nonbipartite theorem is a logical 
consequence of the bipartite one. We describe it here because its proof is substantially easier and 
contains most of the main ideas. Note that Theorem PI does not produce Ramanujan graphs 
because it does not guarantee any control of the least eigenvalue A^. 

We remark that as they are unions of independent matchings, the graphs we produce may have 
multiple edges between two vertices. Thus, they are strictly speaking multigraphs, and do not 
subsume the previous results if one insists on simple graphs. However, it seems that it should be 
more difficult to construct Ramanujan graphs with multiedges than without. Like |MSS15bj . this 
paper establishes existence but does not give a polynomial time construction of Ramanujan graphs. 

1.2 Related Work 

Infinite families of Ramanujan graphs were first shown to exist for m = p + 1, p a prime, in the 
seminal work of Margulis and Lubotzky-Phillips-Sarnak |Mar88HLPS88| . The graphs they produce 
are Cayley graphs and can be constructed very efficiently, and their analysis relies on deep results 
from number theory, which is responsible for the “Ramanujan” nomenclature. Friedman IFriOSj 
showed that a random m—regular graph is almost Ramanujan: specifically, that a union of m 
perfect matchings has non-trivial eigenvalues bounded by 2 y/m — 1 + e with high probability, for 
every e > 0. More recently, in |MSS15b] . we proved the existence of infinite families of m —regular 
bipartite Ramanujan graphs for every m > 3 by proving (part of) a conjecture of Bilu and Linial 
(BL06j regarding the existence of good 2—lifts of regular graphs. Prior to the present paper, it 
was unknown if there exist Ramanujan graphs of every number of vertices. We refer the reader to 
III.WOO and jMSSlhbj for a more detailed discussion of Ramanujan graphs and 2—lifts. 

1.3 Outline of the Paper 

The proofs of both of our theorems follow the same strategy and consist of three steps. In each 
step we present the simpler non-bipartite case first, and then indicate the modifications required 
for the bipartite case. 

First, we show that the expected characteristic polynomials of the random graphs we are inter¬ 
ested in are real rooted and come from interlacing families (reviewed in Section 12.11) , which reduces 
our existence theorems to analyzing the roots of these polynomials. This is achieved in Section [3] 
by decomposing the random permutations used to generate these expected polynomials into swaps 
acting on two vertices at a time, and showing that such swaps correspond to linear transforma¬ 
tions which preserve real-rootedness properties of certain multivariate polynomials. Theorem 13.31 
implies that if A and B are symmetric matrices, then the expected characteristic polynomial of 
A + PBP t is real rooted for a random permutation matrix P. We remark that this argument is 
completely elementary and self-contained, and unlike |MSS15bl IMSS15c| does not appeal to any 
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results from the theory of real stable or hyperbolic polynomials. In the process, we introduce a 
class of “determinant-like” polynomials which may be of independent interest. 

Next, in Section [I] we derive a closed-form formula for the expected characteristic polynomial of 
a sum of randomly permuted regular graphs. We begin by proving that the expected characteristic 
polynomials over random permutations can be replaced by expected characteristic polynomials 
over random orthogonormal matrices. This may be seen as a quadrature or derandomization 
statement, which says that these characteristic polynomials are not able to distinguish between the 
set of permutation matrices and the set of orthogonal matrices; essentially this happens because 
determinants are multilinear, which causes certain restrictions of them to have very low degree 
Fourier coefficients. This component of the proof may also be of independent interest. 

Finally, we appeal to machinery developed in our companion paper (MSSlbaj . which studies the 
structure of expected characteristic polynomials over random orthogonal matrices. In particular, 
such polynomials may be expressed crisply in terms of a simple (and explicitly computable) con¬ 
volution operation on characteristic polynomials, which we call the finite free additive convolution. 
In this framework, the characteristic polynomial of a union of m random matchings is simply the 
m— wise convolution of the characteristic polynomial of a single matching. By applying strong 
bounds on the roots of these convolutions derived in (MS Sinai . we obtain the desired Ramanujan 
bound of 2 \Jm — 1. The requisite material regarding free convolutions is introduced in Sections 12.21 
and 12.31 

These three ingredients are combined in Section [5] to complete the proofs of Theorems 11.11 and 

P 

2 Preliminaries 

2.1 Interlacing Families 

We recall the following theorem from (MSSI5c] , stated here in the slightly different language of 
product distributions. 

Theorem 2.1 (Interlacing Families). Suppose (A;(x)}a;e{o,i} m a family of real-rooted polynomials 

of the same degree n with positive leading coefficient, such that 

-Fji(x) .- ^LJ^flfu (%) 

is real-rooted for every product distribution g = ® ... <g> g m on Q = {0, l} m . Then for every 

k = 1 ,,n and every such g, there is some ojq € Tl such that 

AkC/uJ < Afe(F M ), 

where denotes the kth largest root of a real-rooted polynomial. 

For real rooted polynomials / and g , we write g — > f if the roots of / and g interlace and the 
largest root of / is at least as big as the largest root of g. We will use the following elementary 
facts about interlacing and real-rootedness, which may be found in (Fis08j . 

Lemma 2.2. If g has degree one less than f and both are real-rooted, then 

1. g —> f if and only if f + ag is real-rooted for all a £l 
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2. g —> f implies that f —> f — g. 

If f\ and $2 are monic and real-rooted of the same degree, then they have a common interlacing 
if and only if fi + a /2 is real-rooted for all a > 0. 

2.2 Finite Free Convolutions of Polynomials 

To analyze the expected characteristic polynomials of the random graphs we consider, we will 
need the notion of a finite free convolution of two polynomials, developed in our companion paper 
|MSS15a] , We denote the characteristic polynomial of a matrix by: 

Xx (A) := det(xl - A). 

Definition 2.3 (Symmetric Additive Convolution). Let p(x ) = Xx (A) and q(x) = Xx (B) be two 
real-rooted polynomials, for some symmetric d x d matrices A and B. The symmetric additive 
convolution of p and q is defined as: 

p(x) ffld q(x) = E Xx (-4 + QBQ t ) , 

where the expectation is taken over random orthogonal matrices Q sampled according to the Haar 
measure on 0(d), the group of d-dimensional orthonormal matrices. 

Note that this is a well-defined operation on polynomials because the distribution of the eigen¬ 
values of A -{- QBQ t depends only on the eigenvalues of A and the eigenvalues of B, which are the 
roots of p and q. 

Definition 2.4 (Asymmetric Additive Convolution). Let p(x) = Xx (AA t ) and q(x) = Xx ( BB T ) 
be two real-rooted polynomials with nonnegative roots, for some arbitrary (not necessarily sym¬ 
metric) d x d matrices A and B. The asymmetric additive convolution of p and q is defined as 

p(x) BSE^ q(x) = E X x{(A + QBR t )(A + QBR t ) t ) , 

Q,R 

where Q and R are independent random orthogonal matrices sampled uniformly from 0(d). 

When dealing with a possibly asymmetric d x d matrix M , we will frequently consider the 
dilation 

' 0 M 

M t 0 J ’ 

which is by construction a symmetric 2d x 2d matrix. We will refer to a matrix of this type as 
a bipartite matrix. It is easy to see that its eigenvalues are symmetric about 0 and are equal 
to ±Ai(MM t ) 1/2 , ..., ±A d (MM T ) 1 / 2 , i.e., in absolute value to the singular values of M. This 
correspondence also gives the useful identity 

(mm t ) = Xx ( 

where the operator § is defined by 

(8p)(x) := p(x 2 ). 


0 M 


M 1 


0 


( 1 ) 


5 









With this notation in hand, we can alternately express the asymmetric additive convolution as 


§(p(x) &S d q(x)) = E Xx 

Q,rt 


o 

a t 



Q 

o 


0 

R 


0 

b t 


B 

0 


Q 

o 



( 2 ) 


Explicit, polynomial time computable formulas for the additive convolutions in terms of the 
coefficients of p and q may be found in Theorems 1.1 and 1.3 of [MSS15aj . For this work, we only 
require the following important consequences of these formulas, also established in |MSS15a] . We 
will occasionally drop the subscripts in EB^ and 1++! ^ when it is clear from the context. 

Lemma 2.5 (Properties of EB and I++I ). 1. If p(x) and q(x) are real-rooted then p(x) EB^ q(x) is 

also real-rooted. 


2. If p(x) and q(x) are real-rooted with all roots nonnegative, then p(x) BEEB^ q{x) is also real- 
rooted with all roots nonnegative. 

3. The operations BEb and EEEE^ are bilinear (in the coefficients of the polynomials on which they 
operate) and associative. 

Proof. (1) and (2) are Theorems 1.2 and 1.4 of [MSS 15a] , and bilinearity follows immediately from 
Theorems 1.1 and 1.3 of |MSS15aj . To see associativity, let p(x) = Xx(A),q(x) = Xx(B) and 
r(x ) = Xx (C), and observe that 

(p(i)E?(i))fflr(x) = U^Xx ( QAQ t + RBR t )\ (C) 

\Q R J 

= EE [xx ( QAQ + RBR t ) EB Xx (C)) by bilinearity 
Q R 

= EEExx (QAQ + RBR t + WCW T ) , 

Q RW K ’ 


for random orthogonal matrices Q, R, W. The same argument shows that this is also equal to 
p(x) EB (q(x) EB r(x)). 

An analogous argument using the formula @ shows that I++I is also associative. □ 

A consequence of the above lemma is that for m matrices A\,, A m , identities such as 


® Xx [ 'y ] QiAiQi 

Ql ,..;Qm \ 


Xx (Ai)H3xx (A 2 ) EB - ■ ■ EBXx (A m ) 


(3) 


make sense. 

2.3 Cauchy Transforms 

The device that we use to analyze the roots of finite free convolutions of polynomials is the Cauchy 
Transform. This is the same (up to normalization) as the Stieltjes Transform and the “Barrier 
Function” of |BSSl2l iMSSlhbl IMSSI55] . 
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Definition 2.6 (Cauchy Transform). The Cauchy Transform of a polynomial p(x) with roots 
Ai,..., \d is defined to be the function 


Q P (x) 


i d 

iE 


d x — Xi 
1=1 


1 p'{x) 
d p(x) 


We define the inverse Cauchy Transform of p to be 


/C p (w) = max {a; : Q p (x) = w} . 

Note that the Cauchy transform has poles at the roots of p, and when all the roots A i of p are 
real, Q p (x) is monotone decreasing for x greater than the largest root. Thus, K p (w) is the unique 
value of x that is larger than all the A* for which Q p (x) = w. In particular, it is an upper bound 
on the largest root of p, and approaches the largest root as w —» oo. 

Our bounds on the expected characteristic polynomials of random graphs are a consequence of 
the following two theorems, which are proved in jMSSlhaj . 

Theorem 2.7 (Theorem 1.7 of [MSS15a] ). For real-rooted degree d polynomials p and q and w > 0, 

fc-[El d q M < Kp (w) + £q (w) ~ 1 /W. 

The above theorem is a strengthening of the univariate barrier function argument for character¬ 
istic polynomials introduced in |BSS12| . This may be seen by taking q(x) = Xx (B) = x d ~ 1 (x — d), 
which corresponds to a rank one matrix B = vv T with trace equal to d. It is easy to check that in 
this case p(x) EB q(x) = p(x) — p'{x). 

We remark that Theorem m is inspired by an equality regarding inverse Cauchy transforms of 
limiting spectral distributions of certain random matrix models arising in Free Probability theory; 
we refer the interested reader to |MSS15a] for a more detailed discussion. To analyze the case of 
bipartite random graphs, we will need the corresponding inequality for the asymmetric convolution. 

Theorem 2.8 (Theorem 1.8 of }MSS15 a]). For degree d polynomials p and q having only nonneg¬ 
ative real roots, 

fc§(jttS d q) O) < £§p (w) + /C S q (w) - l/w. 


3 Interlacing for Permutations 

In this section, we show that the expected characteristic polynomials obtained by averaging over 
certain random permutation matrices form interlacing families. The class of random permutations 
which have this property are those that are products of independent random swaps, which we now 
formally define. 

Definition 3.1 (Random Swap). A random swap is a matrix-valued random variable which is 
equal to a transposition of two (fixed) indices s, t with probability a and equal to the identity with 
probability (1 — a), for some a € [0,1]. 

Definition 3.2 (Realizability by Swaps). A matrix-valued random variable P supported on permu¬ 
tation matrices is realizable by swaps if there are random swaps S±, ... ,Sn such that the distribution 
of P is the same as the distribution of the product SnSn-i ■ ■ ■ S 2 S 1 . 
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For example, we show in Lemma 13.51 below that a uniformly random permutation matrix is 
realizable by swaps. 

The main result of this section is that expected characteristic polynomials over products of 
random swaps are always real-rooted. These polynomials play a role analogous to that of mixed 
characteristic polynomials in }MSS15bl IMSS15cj . 

Theorem 3.3. Let A\,... ,A m be symmetric dx d matrices and let {Sjj }i<m,,j<N be independent 
(not necessarily identical) random swaps. Then the expected characteristic polynomial 

Edet ^ n 4 (n j (4) 

is real-rooted. 


An immediate consequence of Theorems 13.31 and 12.11 applied to the family of polynomials 
indexed by all possible values of the swaps Sij, is the following existence result. 

Theorem 3.4 (Interlacing Families for Permutations). Suppose A \,..., A rn are symmetric dx d 
matrices, and P\,..., P rn are independent random permutations realizable by swaps. Then, for 
every k < d: 

< A * (e Xx ^ P % A % Pj 

with nonzero probability. 

Theorem l3.4l is useful because the uniform distribution on permutations and its bipartite version, 
which we use to generate our random graphs, are realizable by swaps. 


A k 


U=1 J 



Lemma 3.5. Let P and S be uniformly random d x d permutation matrices. Both P and P © S 

'P O', 

is the direct sum of P and S. 


are realizable by swaps, where P © S = 


0 5 


Proof. We will establish the claim for P first. We proceed inductively. Let M 2 be a random swap 
which swaps ei and e 2 with probability 1/2, and for k > 2 let 


M k — M k -iSi k M k _i, 


where Si k swaps ei and e k with probability 1/k. 

Let v = (1,2,3,..., d) T . By induction, assume that the first k — 1 coordinates of M k _\v are 
in uniformly random order; in particular, that (. M k -\v)(l ) is a random element of {1,..., k — 1}. 
This means that: 

• With probability 1/k: (. M k _iS\ k M k _\v)(k ) = k and the remaining indices contain a random 
permutation of {1,... , k — 1}. 

• With probability 1 — 1/k: (. M k -iSikM k -iv)(k ) is a uniformly random element j G {1,..., k — 
1} and the remaining indices contain a random permutation of {1 ,... ,k}\ {j}. 

Thus, M k is uniformly random on {1,..., k}, and by induction M c j = P. 

For P © S, we use the above argument to realize P © / and I © S separately and then multiply 
them. □ 





The rest of this section is devoted to proving Theorem 13.31 This is achieved by showing that 
the polynomials in (J4]) are univariate restrictions of certain nice multivariate polynomials. The 
relevant notion is the following. 

Definition 3.6 (Determinant-like Polynomials). A homogeneous polynomial P(X \,..., X m ) of 
degree d in the entries of m symmetric d x d matrices X\.... ,X m is called determinant-like if it 
has the following two properties. 

Hyperbolicity. The univariate restrictions 

q(t) = P(tl - Ai ,..., tl - A m ) 

are real-rooted for all symmetric Ai,..., A m . 

This condition is known as hyperbolicity of the polynomial P(X x,...,A m ) with respect to 
the point (/,/,...,/). We do not discuss the notion of hyperbolicity further, since the self- 
contained definition above suffices for this naDer. We Doint the interested reader to [Pern 12) 
for a detailed discussion of the theory. 

Rank-1 Linearity. For every vector v, index i < rri. and real number s, we have 

P(X 1 ,X 2 ,... ,X{ + svv T ,... ,X m ) = P(X i,... ,X m ) + sD i vv rP{X i,... ,X m ) 

where 

D i,wrP(XX m ) = (J^P(X 1 ,...,X i + svv T , ...,X m )^j | s=0 

is the directional derivative of P in direction (0,... ,vv T ,... ,0), where vv T appears in the 
ith position. Note that D i vv rP(X i,... ,X m ) is homogeneous of degree d — 1. 

An important example of a determinant-like polynomial is the determinant of a sum of matrices: 

P(X 1; ...,X m ) = det(Ai + ... + X m ). 

Hyperbolicity follows from the fact that 

P(tl — Ai ,... ,tl — A m ) = det {mtl — A\ — ... — A m ) 

is the characteristic polynomial of a symmetric matrix. Rank-1 linearity can be seen to follow from 
the invariance of the determinant under change of basis and its linearity with respect to matrix 
entries. Alternatively, one can prove it by using the matrix determinant lemma, which tells us 

det(A'i + svv t + ... + X m ) = det(Ax + ... + A m ) + det(Ax + ... + A m )(Ax + ... + A m ) ■*■). 

The crux of the proof of Theorem 13.31 lies in the fact that random swaps define linear operators 
which preserve the property of being determinant-like. 

Lemma 3.7 (Random swaps preserve determinant-likeness). If P( AT,..., A m ) is determinant-like, 
then for any i < m and random swap S, the polynomial 

K s P(X 1 ,...,SX i S T ,...,X m ) 


is determinant-like. 
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Before proving this lemma, we record some preliminary facts about determinant-like polynomi¬ 
als. 

Lemma 3.8 (Rank-1 updates interlace). Suppose P(Xis determinant-like. Then for 
every vector v and symmetric matrices A\,, A m we have 

P(tl - A 1} ... ,tl - A m ) — P(tl - Ai,... ,tl ~ Ai - vv T , ...tl- A m ), 

where —> denotes interlacing, pointing to the polynomial with the largest root. 

Proof. Assume without loss of generality that i = 1. By rank-1 linearity, 

P(tl - Ai - svv T , ...,tl- A m ) = P(tl - A 1 ,... ,tl - A m ) - sD vvT P(tI - Ai,... ,tl - A m ). 

By the hyperbolicity of P, we know that this is real rooted when viewed as a univariate polynomial 
in t. Since D l vv tP is of degree one less than P, the first part of Lemma 12.21 implies that 

D\ <vv TP(tI — Ai,... ,tl — A m ) —» P(tl — Ai,... ,tl — A m ), 

which in turn by the second part of Lemma 12.21 gives 

P(tl - Ai,...,tl - A m ) —» P(tl - Ai,...,tl - A m ) - D 1 vv tP{A\ -tI,...,A m ~ tl) 

= P(tl - Ai - vv T , ...,tl- Am), 


as desired. 


□ 


Lemma 3.9 (Permutations preserve rank-1 linearity). (1) If Tl is a permutation matrix and 
P(X \,..., X m ) is rank-1 linear then P(LLYiII t , X 2 , ■ ■ ■, X m ) is also rank-1 linear. (2) If P and Q 
are rank-1 linear then so is P + Q. 

Proof. (1) is true because the set of rank one matrices is invariant under conjugation by permuta¬ 
tions. (2) holds because D, vyT is a linear operator. □ 

We will also need the following elementary observation, which says that random swaps corre¬ 
spond to trace zero rank two updates. This is the structural property which causes interlacing to 
occur. 


Lemma 3.10. If a is a transposition and A is symmetric then A — aAa T has rank 2 and trace 0. 

Proof. Assume without loss of generality that a swaps the first two coordinates. Then by symmetry 
the difference A — aAa T has entries 


011 — 022 
021 — 012 
031 — 032 
041 ~ 042 


012 — O21 

022 — On 
032 — O31 
042 — O41 


013 “ O23 
023 — 013 
0 
0 


014 “ «24 
024 — 014 


a P V T 
P —a —v T 
v —v 0 


for some numbers a, /3 and some column vector v of length d — 2. If a / /J then the sum of the 
first two rows is equal to (c, —c, 0,... , 0) for some c ^ 0, and every other row is a scalar multiple 
of this. On the other hand, if a = ft then the first two rows are linearly dependent, and all of the 
other rows are multiples of (1, — 1, 0,... , 0). □ 
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We can now complete the proof of Lemma 13.71 


Proof of Lemma \3.7[ Assume P is determinant-like, and let S' be a random swap, equal to some 
transposition a with probability a and the identity with probability (1 — a). We will show that 

Q(X U ..., X m ) = (1 - a)P(X\,. ..,X m ) + aP(X 1 ,... , aX i( y T , ...,X m ), 

is hyperbolic and rank-1 linear. It is clear that Q(X \,..., X m ) is homogeneous since swaps and con¬ 
vex combinations preserve homogeneity. Lemma [3.91 implies that rank-1 linearity is also preserved, 
so all that remains is hyperbolicity. Assume without loss of generality that i = 1 and consider any 
univariate restriction along (1,1,... ,1): 

Q(tl - A 1 ,... ,tl - A m ) = (1 - a)P(tI - Ai,... ,tl - A m ) + aP(tI - aAia T ,... ,tl - A m ). (5) 

We need to show that this has all real roots. Observe that the second polynomial may be written 
as 

P(tl — Ai — aa T + bb T ,... ,tl — A m ), 

for some vectors a and b , since crAia 1 — A\ is rank two and trace zero by Lemma 13.101 Since P 
is determinant-like, Lemma 13.81 tells us that 

P(tl — A\ + bb T ,... ,tl — A m ) —> P(tl — A\ — aa T + bb T ,... ,tl — A m ) 

and 

P(tl - Ai + bb T , ...,tl- A m ) —■> P(tl - A 1 ,..., tl - A m ), 

whence the two polynomials on the right hand side of (JSJ) have a common interlacing. Lemma 12.21 
then implies that their convex combination must be real-rooted, and the claim is proved. □ 

Applying Lemma 13.71 inductively yields Theorem 13.31 

Proof of Theorem, 1.9. .91 Applying Lemma 13.71 nN times (once for every swap Sjj ) starting with 
P(X i,... ,X m ) = det(]W Xi) tells us that 


n / 1 


N 


E SlJv ...E Sll E 52Jv ...E 5ril det MT [I \ X A 11^ 


i= 1 \j=N 


ki =1 


is determinant-like. Considering the restriction Xi = (t/m)I — Aj finishes the proof. 


□ 


4 Quadrature 

In this section, we show that the expected characteristic polynomials we are interested in are free 
convolutions of the characteristic polynomials of perfect matchings, after the trivial eigenvalues 
corresponding to the all ones vector are removed. This gives us explicit formulas for these polyno¬ 
mials, and more importantly (since we understand the behavior of roots under free convolutions) 
a way of bounding their roots. We begin by showing how to do this for the symmetric case, which 
is more transparent and contains all the main ideas. In Section 14.21 we derive the result for the 
bipartite case as a corollary of the result for the symmetric case. 
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4.1 Quadrature for Symmetric Matrices 

The following theorem gives an explicit formula for the expected characteristic polynomial of the 
sum of two symmetric matrices with constant row sums when the rows and columns of one of 
the matrices is randomly permuted. This can be used to compute the expected characteristic 
polynomial of the Laplacian matrix of the sum of two graphs when one is randomly permuted. In 
this paper, we use the result to compute the expected characteristic polynomial of the adjacency 
matrix when both graphs are regular. 

Theorem 4.1. Suppose A and B are symmetric d x d matrices with A1 = al and 1 = 61. Let 
Xx {A) = (x- a)p(x) and \x (B) = (x - b)q{x). Then, 

EpXx (A + PBP t ) = (x - (a + b))p(x) q{x), (6) 

where P is a uniformly random permutation. 

We begin by writing (JBD in a more concrete form. Observe that all of the matrices A, B,P have 
1 as a left and right eigenvector, which means that there is an orthogonal change of basis V (for 
concreteness, mapping 1 to the standard basis vector e n ) that simultaneously block diagonalizes 
all of them: 

V AV t = A® a, VBV t = B®b, VPV T = P® 1, (7) 

where A® a denotes the direct sum 

'A O' 

0 a 

Since the determinant is invariant under change of basis, we may write 

E P det(xJ -A- PBP t ) = E P det(sJ - VAV T - {VPV T ){V BV T ){V P T V T )) 

= Epdet(x/ — (A © a) — (P © 1 )(B 0 b)(P T 0 1)) 

= (x — a — 6)Ep det(xl — A — PBP t ). (8) 

Notice also that p(x) = Xx (A^J and q(x) = Xx (^B^j , so 

p(x) EB q(x) = Eg det (xl — A — QBQ T ), 

where Q is a (Haar) random {d— 1) x {d— 1) orthogonal matrix. Thus, (j6]) is equivalent to showing 
that 

E^ det {xl -A- PBP t ) = Eq det(xl -A- QBQ T ), (9) 

for all {d— 1) x {d— 1) symmetric matrices A, B. Note that for any permutation P, the orthogonal 
transformation P correspondingly permutes ej,..., e n , the projections orthogonal to 1 of the stan¬ 
dard basis vectors ei,..., e^, embedded in M d_1 . Since these are the vertices of a regular simplex 
with d vertices in centered at the origin, we interpret the P as elements of the symmetry 

group of this simplex. We denote this subgroup of 0(d — 1) by A^-i- 

Since there is no longer any assumption on A, B other than symmetry, we may absorb the xl 
term into A in (f9j) . and we see that it is sufficient to establish the following. 
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Theorem 4.2 (Quadrature Theorem). For symmetric d x d matrices A and B 


E det (A + PBP t )= E det (A + QBQ T ). 

P&A d QeO(d) 


( 10 ) 


It is easy to see that the theorem will follow if we can show that the left hand side of (HOD is 
invariant under right multiplication of P by orthogonal matrices. 

Lemma 4.3 (Invariance Implies Quadrature). Let f be a function from 0{d ) to M and let H be a 
finite subgroup ofO(d). If 

( 11 ) 


E f(P) = E f(PQ o). 

P&H v ’ P£H v ' 


for all Qo € O(d), then 


E f(P)= E f(Q), 

P&H J Q&0(d) 


( 12 ) 


where Q is chosen according to Haar measure and P is uniform on H. 
Proof. 


E f(Q)= E E f(PQ)= E E f(PQ) 
QeO(d) QeO(d) P&H p&h QeO(d) 

= E E f(P) = E f(P), 

P£H QeO{d ) P&H 


as desired. □ 

We will prove Theorem 14.21 bv showing that f(P) = det {A + PBP 1 ) satisfies (HU). We will 
achieve this by demonstrating that / is invariant under certain elementary orthogonal transforma¬ 
tions acting on 3-faces of the regular simplex, which generate all orthogonal transformations. Let 
us fix some notation to precisely describe these elementary transformations. 

Given three vertices ii,ej,ek of the regular simplex, let Aij^ denote the subgroup of Ad con¬ 
sisting of permutations of ei,ej,ek which leave all of the other vertices fixed. Let Oi jk denote 
the subgroup of 0(d) acting on the two dimensional linear subspace parallel to the affine subspace 
through the three vertices, and leaving the orthogonal subspace fixed. Note that A % ] k is a subgroup 
of Oi jk, an d that these groups are isomorphic to A 2 and 0(2), respectively. 

The heart of the proof lies in the following lemma, which implies by Lemma 14.31 that the 
polynomials we are interested in are not able to distinguish between the uniform distributions on 
A 2 and 0(2). The reason for this is that these polynomials have very low degree (at most two) 
in the entries of any orthogonal matrix Q acting on a two-dimensional subspace, a fact which is 
essentially a consequence of the multilinearity of the determinant. The argument below is similar 
to the proof of Lemma 2.7 in |MSS15a] . 

Lemma 4.4 (Invariance for A 2 ). If A and B are symmetric d x d matrices, then for every Q 0 € 

0 ( 2 ), 

E det (4 + (P © I d - 2 )B(P © I d _ 2 ) T ) = E det (A + ( PQ 0 0 I d - 2 )B(PQ 0 © I d - 2 f). (13) 
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Proof. Let 50(2) be the subgroup of 0(2) consisting of rotation matrices 


Rg — 


cos tt sm tt 
— sin 6 cos 9 


and let Z 3 be the subgroup of A 2 consisting of the three rotations R t ,t G T := {0, 27t/3, 47t/ 3}. 
We begin by showing that 

E det(A + (P®I)B(P®I) T ) = E det (A + (PRg ® I)B(PR e ® if), (14) 

P&Z3 PGZ3 

for every 9, where / is the (d — 2)-dimensional identity. Since the elements of Z 3 are themselves 
rotations, we can rewrite thrice the right hand side of (11411 as 

det (A + ( R T Rg © I)B{R t Rq © if ) 

r£T 

= y ^ det (A + (R r +g ffi I)B{R r+ e © if) 

reT 

2 

= for some coefficients Ck, by Lemma 1431 

tGT k =—2 
2 


c k e ike (e m + e ik27r / 3 + e ifc47r / 3 ) 


k =-2 

= 3co since the terms with \k\ = 1,2 vanish. 

As this quantity is independent of 9, we can assume 9 = 0, which gives the left hand side of (1141) . 
To finish the proof, we observe that 

E det(A + (P ® I d _ 2 )B(P ® I d _^) T ) = E E det (A + {PD ffi I)B(PD ffi I) T ) 
pga 2 d&fp&z 3 

= E E det (A + (P ffi 7)(D ffi I)B{D ffi I) T (P ffi I) T ), 

DeF P£Z 3 

where F consists of the identity and the refiection across the horizontal axis: 


F := 


1 0 
0 1 


1 0 
0 -1 


and D is chosen uniformly from F . 

Thus, the left hand side of (|13l) is invariant under conjugation of B with the matrices D®I,D G 
F. Since every Qq G 0(2) can be written a RgD for some D G F, and we have already established 
invariance under Rg ffi I in (1141) . the lemma is proved. 

□ 


Lemma 4.5 (Determinants are Low Degree in Rank 2 Rotations). Let A, B be d x d matrices and 
let 

det (A + (R e ffi I d - 2 )B(R e ffi I d _ 2 f ) = Y c k e iM . 

k 

Then c k = 0 for \k\ > 3. 
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Proof. Recall that all 2 x 2 rotations may be diagonalized as 


Re 


cos 9 

sin 9 

= U 

e 16 

0 

— sinf? 

cos 9 

0 

e -ie 


where 

1 1 ' 

i —i 

is independent of 9. This implies that ( Rq (B Id- 2 ) = VDV^ for diagonal D containing e lB and e~ lB 
in the upper right 2x2 block and ones elsewhere, with V independent of 9. Thus, we see that 



det (A + (Rq 0 Id-2)B(Re 0 Id- 2 )^) — det (A(Re 0 Id-2 ) + (Re 0 Id-2)B) 

= det ( AVDV t + VDV^B^j 

= det (v^AVD + DV^Bv'j 

Notice that the matrix M = V^AVD + DV^BV depends linearly on e™ , and that the e lB 
(resp. e~ lB ) terms appear only in the first (resp. second) row and column of M, respectively. Since 
each monomial in the expansion of the determinant contains at most one entry from each row and 
each column and e lB ■ e~ lB = 1 , this implies that no terms of degree higher than two in e lB or e~ lB 
appear. □ 


Corollary 4.6 (Invariance for Aij^)- For every i, j and k, 

E det (A + PBP t )= E det (A + QBQ T ). 

PeA i:jjk Q&O i:j , k 

Proof. Let V be the orthogonal transformation that maps the affine subspace spanned by the 
vertices e t . e.j , to the first two coordinates of M 2 , with any one vertex mapped to a multiple of e±. 
Conjugation by V maps Ajj ^ to A 2 0 Id -2 and Oi.j.k to 0(2) 0 Id- 2 , abusing notation slightly in 
the natural way. Since the determinant is invariant under change of basis, Lemma 14.41 tells us that 

E det (A + PBP t ) = E &et(V AV T + (P 2 0 I)VBV T (P9 0 I) T ) 

P&A iijik P 2 &A 2 

= E det (VAV t + (Q 2 ®I)VBV t (Q 2 ®I) T ) 

Q2&0(2) 

= E det (^4 + QBQ t ), 

Q£-Oi,j,k 

as desired. □ 


Lemma 4.7 (Ojj ^ generate 0(d)). Given a regular simplex in M. d , the union over i, j, and k of 
Oi,j,k generates 0(d). In particular, every matrix in 0(d ) may be written as a prodiLct of a finite 
number of these matrices. 

Proof. Let T^ be the subgroup of 0(d) generated by Uj j k Oi,j,k<h- Let eo,..., be the vertices 
of the regular simplex. For 1 < h < d, let Eh be the linear subspace parallel to the affine subspace 
through eo,... ,ih- We will prove by induction on h that T h contains the action of the orthogonal 
group on Eh . The base case is h = 2, for which Oo,i ,2 is precisely the action of the orthogonal 
group on E 2 . 
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Assuming that we have proved this result for h — 1, we now prove it for h. To this end, let 
Uh = &hi and let u\,..., Uh-i be arbitrary orthonormal vectors in Eh that are orthogonal to Uh■ We 
will prove that for every orthonormal basis w \,... ,Wh of Eh, there is a Q £ Th such that Qwi = U{ 
for 1 < i < h. 

We first consider the case in which w^eh > 0. Let Fh denote the 2-dimensional affine subspace 
spanned by {eh, eh-i,ih- 2 }, and observe that there must be a unit vector p € EhH Fh with p T ih = 
wj{eh- This follows because the intersection of Fh with the unit sphere in Eh is a circle containing 
{eh, eh- 1 , eh- 2 }, P H > p T eh is a continuous function, and we have ej{eh = 1 and e^^e-h = e%_ 2 eh < 0. 
As eh is orthogonal to Eh -1 and ih is invariant under T/ 1 _ 1 , the induction hypothesis implies that 
there must be a T € Th_i so that Twh = p■ Moreover, there is an element T 2 of Oh- 2 ,h-\,h that 
maps p to eh- So, their composition W = T 2 T sends Wh to eh- Since W is orthogonal, it must 
send w\,... ,Wh -1 to Eh- 1 , and so by induction may be composed with a map in Th-i that sends 
Ww\,..., Wwh-i to tti,..., Uh -1 without moving eh- The resulting map is the desired Q. 

In the case that w^eh < 0, we begin by applying a map in Th that sends Wh to a vector that is 
orthogonal to eh so that we can then apply the previous argument. For example, we can do this 
by defining p to be one of the two unit vectors in Fh with p T e-h = —wj{eh- We then apply a map 
in T/j_i that sends Wh to —p, and then a map in Oh- 2 ,h-i,h that maps p, and thus also —p, to a 
vector orthogonal to eh- □ 


Theorem 4.8 (Invariance for Aj). Let A and B be d x d matrices, and let 

Ia,b{Q) = det (A + QBQ t ) . 

Then, for all Qo € 0(d), 

E f AB (P) = E }a,b{PQo)- 

P&A d P&A d 

Proof. We will use the fact that 

E f A>B (P) = E E f A , B {PP2) = E E f P T APtB (P 2 ). 

PeA n PeA d p 2 eA itjtk P&A d Pi&A itjtk 

Applying Corollary 14.61 reveals that for every Q 2 € 

E. E W, S (P 2 ). = E E fp T AP, B (P2Q2)- 

P&A d P 2 &A i: j :k P&A d P 2 &Aij ik 

= E E f A , B {PP2Q2). 

P£A d P 2 eAi tjtk 


Thus, we conclude that 


= E f A ,B(PQ2)- 

P&A d 


E f AiB {P) = pE f A , B (PQ2) 

P&A d P&A d 


for every Q 2 € O itjik , for every i,j, k. 

Let Q 0 € O(d). By Lemma 14.71 there is a sequence of matrices Q 1 ,..., Q m , each of which is in 
Oi,j,k for some i, j and k, so that 

Q 0 — Q1Q2 ' ' ' Qm- 

By applying the previous equality m times, we obtain 


E f(PQo) = E f(PQi ■ ■ ■ Q m ) = E f(P). 

P&A d PeA d PeA d 


□ 
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Proof of Theorem \ f.2\ Follows from Theorem 14.81 and Lemma 14.31 □ 

Proof of Theorem \f.l\ Follows from Theorem 14.21 Q , and ([ 8 ]) . □ 

We conclude the section by recording the obvious extension of Theorem 14.11 to sums of m 
matrices. 

Corollary 4.9. Let A\,... ,A m be symmetric d x d matrices with A,1 = afl and Xx (Ai) = (x — 
a-i)pi(x). Then, 

PiAP^j = 

where P \...., P m are independent uniformly random permutation matrices. 

Proof. Apply a change of basis so that each A* = A % © a*, divide out the (x — YllLi a *) term as in 
([8]) . and apply Theorem 14.21 inductively (to — 1) times, replacing each Pi with a random orthogonal 
Qi (this requires conditioning on the other Pj and Qj , but by independence the distribution of each 
Pi is still uniform on Af). Finally, appeal to the identity Q to write this as an m— wise additive 
convolution. □ 


m \ 

X ~ ^2 ai ) Pl(x)S...Spm(x), (15) 


E 

Pl,-,Prr 


Xx 


E 

\i =1 


4.2 Quadrature for Bipartite Matrices 

Theorem 4.10. Suppose A and B are (not necessarily symmetric) d x d matrices such that A1 = 
A t 1 = al and B1 = B 1 1 = bl. Let Xx {AA T ) = (x — a 2 )p(x) and Xx {PP T ) = (x — b 2 )q(x). 
Then, 

p E T'( 

= (x 2 - (a + 6 ) 2 )S (p(x) BEEEq(x)) , (17) 


0 

a t 


A 

0 


+ {P 0 S) 


0 

b t 


B 

0 


(P © S) T j = § ((x - (a + b) 2 )p{x) EEHg(x)) 


(16) 


where P and S are independent uniform random permutation matrices. 

As in the nonbipartite case, we begin by applying a change of basis V that isolates the common 
all ones eigenvector and block diagonalizes our matrices as: 

VAV t = A © a, VBV t = B © b, VPV T = P © 1, VSV T = S © 1. (18) 


Conjugating the left hand side of (fl 6 l) by (V © V), we see that it is the same as 
0 


E Yjl J .rp 

p,s VL(^©o) t 


= E Xx 

P,S AX 


(.A © a) 

0 

0 


+ ((P© 1) © (S'© 1)) 


0 (P © b) 
(B © b) T 0 


_{A + PBS t © (a + b)f 

T 


(A + PBS t © (a + b)) 
0 


= E § X x ( (A + PBS 1 © (a + b))(A + PBS 1 © (a + b)) J 
P, s 


T\T 


= (x 2 - (a + by) E S Xx ({A + PBS 1 )(A + PBS 1 ) 

PjS 


= (x 2 -(a + b) 2 )E s Xx 


' 0 A 

A T o 


+ (P©S) 


0 B 
B t 0 


(P©S ) 5 


((P©1)©(S©1)) J 


(19) 
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As in the previous section, the matrices P and S are random elements of the group Ad- 1 - Observe 
that 

p{x) = Xx (AA T ) , and q(x) = Xx(BB T ' 


Recalling from <[2|) that 

§ (p(x) @EB q(x)) = E QijRe 0 (d-i)Xx 


0 A 
A t 0 


+ (Q © R) 


0 B 
B t 0 


(Q 0 R) J 


and removing all the ~s as before to ease notation, we see that the conclusion (1161) of Theorem 14. 101 
is implied by the following more general quadrature statement. 


Theorem 4.11. For all symmetric 2d x 2d matrices C and D: 

E Xx(C + (P®S)D(P®S) t )= E Xx(C + (Q ® R)D{Q ® R) t ) 

P,SeA d Q,RGO(d) 


( 20 ) 


This theorem is an immediate consequence of two applications of the following corollary of 
Theorem 14.21 from the previous section. 

Corollary 4.12. If C and D are symmetric 2d x 2d matrices, 

E det(C + (P © I)D(P © I) T ) = E det(C + (Q © I)D(Q © I) T ). 

P&A d QeO(d) 

Proof. The proof is identical to the proof of Theorem 14.21 except we replace P £ Ad with P © I 
and Q £ 0(d) with Q © I at each step. 

Specifically, let 

fcAQ) ■= det {C + {Q © I)D{Q © I) T ). 

Applying Corollary 14.61 as before reveals that for every i,j, k and every Q 2 € 

E f c , D (P) = E p E fcAPPz) = E p E fcAPPM = P E fcAPQ*)- 

P^Ad P^Ad P^zAd P2GAij t k P^Ad 

Since an arbitrary Q 0 € 0(d) is a finite product of such Q 2 by Lemma 14.71 this means that 


p E d f cM PQ B ) = E A J c MP) 

for all Q 0 G O(d). Lemma l4~3l completes the proof. 

Proof of Theorem \f.ll\ Since P and S are independent, we have 
E XxiC + {P®S)D{P®S) t ) 

P,seA d v ’ 

= E E det(xI + C + {P ©/)(/© S)D(I®S) t (P®I) t ) 

SGA d PGA d 

= E E det (xI + C+(Q®I)(I®S)D(I®S) t (Q®I) t ) by Corollary [442] 
SeA d QgO(d) 

= E E det (xl+ (Q®I) T C(Q®I) + (I®S)D(I®S) T ) 

Q^O(d) S^Ad 


□ 


E det(sJ + (Q © lfC{Q ©/) + (/© R)D{I © R) T ) by Corollary WM 


= E 

QeO(d) R&0(d) 

= E det (xI + C + (Q®R)D(Q®R) T ), 

Q,R&0{d) 


as desired. 


□ 
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Proof of Theorem \4-10\ Follows from Theorem 14.111 (fl8l) , and (fl9l) . 

As before, Theorem 4.10 extends effortlessly to the case of many matrices. 


□ 


Corollary 4.13. If A\,... ,A m are matrices with Ai 1 = Aj 1 = ai and Xx {AiAj) = (x-af)pi(x), 
then 


E Xx 


0 Ai 

AT 0 


(Pi © Si f 


= I x - 


\i=1 

y , j j § \pi(x) E3 ... E3p m (x)], 


where the Pi and Si are independent uniformly random permutations. 
We omit the proof, which is identical to the proof of Corollary 14.91 


5 Ramanujan Graphs 


In this section, we combine the Cauchy transform, interlacing, and quadrature results of the previous 
sections to establish our main Theorems 11.11 and P 

Proof of Theorem \1.2[ Let M be the adjacency matrix of a fixed perfect matching on d vertices, 
with d even. Since the uniform distribution on permutations is realizable by swaps (Lemma 13.5[) . 
Theorem 13.41 tells us that with nonzero probability: 


A 2 PiMPj < A 2 E Xx £ PiMP} 


K i— 1 


\i— 1 


Corollary 14.91 reveals that the polynomial in the right-hand expression may be written as an 
m-wise symmetric additive convolutiord 


E(x) := E xx p i A i p I = ( x - m)p(x) l+ld-i • • • SU-i p(x) (m times), 

•* !}•••>-* TTl 


\i= 1 


where 

P (x)== ( x _ i) d / 2 ~\ x +1)^ 2 , 
x — 1 

is the characteristic polynomial of a single matching with the trivial root at 1 removed. Our goal 
is therefore to bound the largest root of p(x) EE... EE p(x), which is the second largest root of E(x). 

We will do this using the inverse Cauchy transform described in Section 12.31 The Cauchy 
transform of p(x) is given by 


[x = 


d/2-1 1 


d/2 1 


d — 1 x — 1 d — 1 x + 1 


2 We remark that this formula works for arbitrary regular adjacency matrices Ai, as mentioned in the abstract. 
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Notice that for every x > 1, putting the trivial root at 1 back only increases the Cauchy transform: 


G P (x) < 


d/2 1 d/2 1 

d x — 1 d x + 1 


x 2 _ i g x(M) ( x ) • 


( 21 ) 


Since both functions are decreasing for x > 1, this implies that the inverse Cauchy transform of p 
is upper bounded by that of %(M): 


JC P (w) < /C X ( M ) (w ), 


for every w > 0. 

Applying the convolution inequality in Theorem 12.71 (m — 1) times yields the following upper 
bound on the inverse Cauchy transform of the m —wise convolution of interest. 


{w)<m- K p (w) - < m • /C x(M) (w) ~ m 


w 


w 


( 22 ) 


Recalling from (121|) that 


x 


K x m (w) = x <<=► w = x2 _ l . 
the right hand side of (|22D may be written as 


771—1 (t 77 — l)(x 2 — 1) X 2 + (?77 — 1) 

mx -= mi-=-. 


w 


x 


x 


which is easily seen to be minimized at x = y/m — 1 with value 2 y/m — 1. Thus, the second largest 
root of E{x) is at most 2y/m — 1. □ 


Proof of Theorem 11.11 Let 


M = 


0 I 
I T 0 


be the adjacency matrix of a perfect matching on 2d vertices, across the natural bipartition. Then, 
for independent uniformly random dx d permutation matrices Pi,..., P m , Si ,..., S m , the random 
matrix 

A = j2(p i ®s i )M(p i ®s l ) T = jr (p f jT) 

is the adjacency matrix of a union of m random matchings across the same bipartition. Since the 
distribution of the (Pj © Sf) is realizable by swaps iLemma 13.51) . Theorem 13.41 implies that 


HA) < A 2 ( E X x ( 0 HM(P, 0 Si) 


\ i— 1 


with nonzero probability. Since 1 1 = 1, Corollary 14.131 implies that the polynomial on the right 
hand side is equal to 


(x 2 — 772, 2 )S[p(a:) EEE rf-i ... I H—H d_ 1 p(x)} (m times), 


where 


V{x) = Xx (Id-ilj-i) =( x ~ 1 ) 


d -1 
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We upper bound the inverse Cauchy transform of this m— wise convolution using Theorem 12.81 


fc§(pE3...E3p) ( w ) < m ■ ICsp (w) 


m — 1 
w 


m • K. {x i_ x)d -1 ( w) 


m — 1 
w 


Since 

X 

^{ x 2 — l) d_1 (^0 “^2 Y ’ 

this is now identical to the calculation (1221) . so we obtain again the bound 2 y/m — 1. Thus, we 
conclude that A 2 CA) < 2-\Arc — 1 with nonzero probability. Since A is bipartite, its spectrum is 
symmetric about zero, so we must also have Ad_i(A) > —2 \Jm — 1, whence A is Ramanujan. □ 
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